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Abstract 

Operating just once the naive Foldy-Wouthuysen-Tani transformation on 
the Schrodinger equation for Qq bound states described by a hamiltonian, 
we systematically develop a perturbation theory in l/mq which enables one 
to solve the Schrodinger equation to obtain masses and wave functions of the 
bound states in any order of l/mq. It is shown that positive energy projection 
with respect to the heavy quark sector of a wave function is, at each order 
of perturbation, proportional to the 0-th order solution. There appear also 
negative components of the wave function except for the 0-th order, which 
contribute also to higher order corrections to masses. 
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I. INTRODUCTION 

Since the introduction of the ingenious notion of heavy quark effective theory 
many physical quantities, especially regarding to B mesons, have been calculated. Further- 
more formulation has also been polished to incorporate higher orders.! However, the way 
to incorporate higher order terms is somewhat insufficient since only the operator prod- 
uct forms in higher orders are obtained but their matrix elements are just parameters and 
must be determined so that the physical quantities are fitted with experiments.! Another 
obscure point of HQET is that the residual momentum, which is a difference between a 
meson momentum and a heavy quark momentum, necessarily comes in and we do not know 
how properly this quantity can be taken into account as a perturbation. On this point 
HQET is different from QED in the sense that the next order perturbation terms cannot be 
determined automatically from the lower orders. 

Let us remind how HQET has been constructed to calculate physical quantities regard- 
ing heavy mesons/baryons. When HQET was originally developed, lagrangian formulation 
was adopted and light/heavy quarks are treated as elementary fields. In this approach an 
appropriate transformation only on heavy quarks is introduced so that heavy quark fields 
are separated into positive /negative components in energy. This picture was intuitively be- 
lieved to realize the situation that a positive component of a heavy quark behaves like a 
static color source and a light anti-quark with gluon clouds is relativistically hanging around 
it, which owes most of the messy part of the strong interactions. To stress even more this 
intuitive idea, one assumes that a heavy quark carries most of the momentum of a heavy 
meson, which yields a residual momentum. 

Even though the above formulation has been well developed, this cannot be applied 
straightforwardly to the Schrodinger equation for a bound state of a Qq system. It is 
the purpose of this paper to construct a perturbation theory in l/mg a la HQET for a 
hamiltonian system, i.e., for a bound system which is described by a hamiltonian with 
an appropriate potential. Instead of particle fields, wave functions of a Qq system at rest 
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are the objects to deal with. Then matrix elements of currents/operators between heavy 
mesons are calculable at any order in principle since meson wave functions are calculated 
as solutions to the Schrodinger equation together with their masses as eigenvalues with 
corrections and there appears no residual momentum of a heavy quark in this formulation. 
In our formulation, at first the naive Foldy- Wouthuysen-Tani (FWT)@ transformation is 
introduced and instead of working only with positive components of heavy quark fields there 
appear negative components as well and hence the formulae to calculate those negative 
components are also given at each order of perturbation. This formulation gives diffrent 
results from the previous calculations! on the masses of D , D* , B and B* because of negative 
components of heavy quark sector which is briefly mentioned at the end of this letter. 



II. HAMILTONIAN 

A projection operator which extracts a positive component of a heavy quark plays an 
important role in HQET. In a much more precise statement this projection operator is 
defined so that the equation for the two upper components of a Dirac spinor (positive 
components) does not mix with the two lower components (negative components) in the 
lagrangian. Here before introducing a projection operator, let us recall what the naive FWT 
transformation does. This transformation, when it operates on a free Dirac Schrodinger 
equation, turns a kinematic hamiltonian, a ■ p + mj3, into EjS and leaves us only the two 
upper components of a free wave function. 

UpwT {c( ■ p + m(3) Up^T Ufwt4> = E Ufwt<P E (3 Ufwt4> = E Ufwt4>, 

i.e., 

E{l-P) Ufwt(P = 0, 

which means that (1 + P)/2 UFWT(f> = Ufwt4>^ only the two upper components of 
Ufwt4> survive. Here 
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U FWT (p) — exp {W (p) ^ ■n) — cos W -\- ^ ■ n sinW, (1) 

As one can see, in this formulation a projection operator, (1 — /3)/2, naturally comes in to 
exclude two lower components. Based on this fact, introduction of the FWT transformation 
may give us a good starting point for studying a Qq system as in fact will be seen later 
rather than trying to develop a perturbation theory of the original hamiltonian without any 
t r ansf ormation . 

The state of a bound state in this paper is defined by 

IV') - jd'^j d'y ^ap{x - y) ql \x) Q\{y) |0) , (3) 

where is a charge conjugate field of a light quark q XSand its conjugate state by {tp\ = 
with (0| = |0)^ Given a hamiltonian density, 7i, for a Qq system, varying the quantity 

{i;\{n-E)\i;), 

in terms of ipapix — y), and setting it equal to zero, then we have the matrix form of the 
Schrodinger equation 

H i/j = (mq + E) ^. 

Here a bound state mass E is broken into two parts, tjiq and the rest, E. A binding energy 
is given by E — rUq. Operating the FWT transformation only on a heavy quark sector in this 
equation at the center of the mass system of a bound state, one can modify the Schrodinger 
equation as, 

{Hfwt - mq) ipFWT = E i/jfwt, (4) 

where 

Hfwt = Ufwt (p'q) HUp^^ (pq) > i^FWT = Ufwt [pq] ip- (5) 

Note that the argument of the FWT transformation operated on a hamiltonian from left is 
different from the right-operated one, since an outgoing momentum, pq ', are different from 



incoming one, pq. Expressing model-dependent instantaneous potential terms as iJ'"*, a 
hamiltonian is given by 

H ^aq-pq + mqf3q + ■ PQ + UIq/Jq + . (6) 

Next HpwT — TTT'Q can be expanded in I/ttiq, which is given by 

Hfwt -rnq^ H^i + Ho + Hi + H2 + ---, (7) 



where 



H_^ = rriQ {I3q - 1) , (8a) 

Ho^aq- p + mq(3q + Hp^To, (8b) 

Hr = -^p'f3Q + W;^^^„ (8c) 

H2 = -f^F;^r2) (8d) 



Here stands for the z-th order hamiltonian, H'^p^n ^'^^ the i-th order of the interaction 
terms and since a bound state is at rest, 

P = Pq = -PQ, p'=Pq' = -PQ', Q=P'-P, 
are defined, where primed quantities are final momenta. 

III. PERTURBATION 

Using the hamiltonian obtained in the last section, we show in this section that the 
Schrodinger equation can be solved order by order in l/mq utilizing a projection operator. 
First we introduce projection operators: 

A. ^ ii*. (9) 

which correspond to positive-/negative-energy projection operators for a heavy quark sector 
at the rest frame of a bound state. These are given by (1 ± ^)/2 in the moving frame of a 



bound state with = P^/E in HQET where P'^ is the four-momentum of a bound state. 
Then we expand the mass and wave function of a bound state in l/mg as 

E = E^ + E\ + E^ + ..., (10) 
i^FWT = ipo + V? + + • • • , (11) 

where k stands for a set of quantum numbers that distinguish independent eigenf unctions 
of the lowest order Schrodinger equation, and a subscript i of E^ and for the order of 

A. -1st order 

From eqs.(|D and (pal), the -1st order Schrodinger equation in l/mg is given by 

- 2mQA_^o' = 0, (12) 

which means 

^0 = A+^o'- (13) 

That is, the 0-th order wave function has only a positive component of the heavy quark 
sector whose form is given below. 

B. 0-th order 

The 0-th order equation is given by 

Multiplying projection operators, A±, from left, respectively, we obtain 

-2mQA_V'i' + A_/Jo^o = 0. 
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(14) 

(15) 
(16) 



Eq. (|T5D gives the lowest non-trivial Schrodinger equation. Let us assume this equation 
is solvedl and all the independent eigenfunctions and eigenvalues are obtained analytically 
and/or numerically as and 

/ufc(r,n) 0\ 

V^o = -^-^ = , (17) 

where Uk and are 2 by 2 matrices. Likewise we can in general assume that a A_ component 
of any wave function is expanded in terms of 

^k = [ , 1 ■ (18) 

Expanding in terms of this set of eigenfunctions, one can solve Eq. ([T6| ) as follows. 

Setting 

A_^/.f = ^cf$7, (19) 
I 

one obtains coefficients, c", as 

= / r'^drd^ tr k.E^k^^f) , (20) 

where the 0-th order wave function is normahzed to be 1, 

Jd'rtT{^f^f)=5,,. (21) 

C. 1st order 

The 1st order equation is given by 

- 2mQA_V^2^ + ifo^f + = Eo^^^ + E'l^. (22) 

As in the above case, multiplying projection operators, we obtain 

A+Hoiji + A+i^iV-o = ^0 A+^i' + i^iVo , (23) 
-2mQA_V'2 + A^Ho^i + A-ifiV-o = A_^^ (24) 
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When one looks at Eq. (^), one easily notices that this can be separated into two indepen- 
dent equations. One equation gives a solution to the positive component of ip^ as 

which means from Eq. (|T3|) that 

A+^^ = ct+$+, (25) 

and the coefficient, cf_,_, should be determined by a normalization of the wave function up 
to the first order, 

c^V tr(^/'^V^) = 1, (26) 

whose definition is allowed because here we are not calculating the absolute value of the 
form factors. The appropriate normalization will be determined in future papers in which 
we will give several kinds of form factors. Using this equation, one easily notices that up to 
the first order, 

= 0, (27) 

where the coefficient is assumed to be real. This completes the solution for tpf since A^ipf 
is obtained in the last chapter, i.e., 

i^i=T.c{'-<^7- (28) 

e 

Another equation is given by 

A+/JoA-^i' + A+H,4 = E'[4, (29) 

which gives the first order perturbation corrections to the mass when one calculates matrix 
elements of the Ihs among eigenf unctions, $^ like in the ordinary perturbation of quantum 
mechanics. 

Eq. ( P^ ) gives a A_ component of 1^2 case of A^ipf in subsection [illB| , i.e., 

setting 
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(30) 



one obtains coefficients, 



Jk 

-2- 



2mQ 



r drdQ tr 



(31) 



D. 2nd order 

The 2nd order equation is given by 

- 2mQA_^lJ^ + H^^l + H^^'l + i^a^o = ^oVs + ^iV? + ^2^0- (32) 
As in the above cases, multiplying projection operators, we obtain 

A+i7o^2 + A+H,i,\ + A+H^i^l = Eo'A+V^2 + ^1 A+^f + ^a'^o, (33) 
-2mQA_V^^ + A_i/oV^2 + A^i^iV^f + A./^s^o = ^0 A_^2^ + A.^^f . (34) 

When one looks at Eq. (|33D, one easily notices that this can be separated into three inde- 
pendent equations. The ffist one gives a solution to the positive component of ^/'f : 

A+HoA+^Pl = E^A+^l 

which means from Eq. (|T5[) that 

A+^2' = (35) 

The coefficient, C2+, should be determined by a normalization of the wave function up to 
the second order using Eq. (pU]) and by assuming the coeffcient is real, it is given by 




There are some other contributions to this coefficient, C2+, from transitions of $^ to other 
states and back to $^ which are obtained when one diagonalizes the whole energy matrix. 
This completes the solution for ■j/'f when combined with the A_ component obtained above. 
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=ct'f^ + E4-'f7- (37) 

I 

The second equation is given by 

aA+i7iA_^t' + A+^iA+^i = A+V-i , 

where a is determined so that this equation coincides with the first order energy correction 
equation, Eq. (pOf), using A+^z^^ = c^^i/jq-. 

a = c'l^ = 0. 

The third equation is given by the remainder, 

A+HoA^i;!^ + A+i/iA_^t + A+i^sV^o = E^Ki^t (38) 

which gives the second order corrections to the energy by taking the matrix elements of the 
Ihs among 

Eq. ( P^ ) gives a A_ component of 1IJ3 as in the cases of A_iljf and A.-j/'l given in subsec- 
tions [111B| and |ill C[ i.e., setting 

A-^HE^-'^Z, (39) 

e 

one obtains coefficients, 

$7^ [A^Hoij^ + A^Hiij^ + A_i/2$fc - ^0 A_^2^ - A„V^f)" . 

(40) 



r drdVt tr 



E. Higher order 

The i-th order equation is given by 

- 2mQA_^t, + ^ if,^^. = E^^lj- (41) 

j=0 j=0 

We follow the previous subsections to derive all the results from this equation below. Mul- 
tiplying projection operators, we obtain 
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EA+i/,#_,. = E4^+^'-., (42) 

j=0 j=0 

-2mQA_V^ti + E ^~H,4jI^ = (^3) 

j=0 j=0 

Eq. ( ^21) gives a solution to a positive component of ■i/'f as 

A+i/oA+^- = i^oA+V'f , (44) 

which means 

A+^f = ct$^, (45) 

where the coefficient, c*^_, should be determined by a normalization of the wave function, 
Eq. up to the i-th order and there are of course other contributions mentioned in 

subsection |III D| . This completes the solution for when combined with the A_ component 



obtained from the {i — l)-th order equation. Eq. (|4^) also includes other equations, which 
are identified as the {i — j)-th order (j = 1 ~ i — 1) energy correction equations. These 
are straightforwardly obtained order by order though complicated and hence are omitted 
here. After extracting these equations, the remainder of Eq. (^) gives the z-th order energy 
correction equation. Eq. ( ^31) gives a A_ component of ipf^i, i.e., setting 

A-^ti = Ec"i,-'^7, (46) 

i 

one obtains coefficients, c";^ _, as 



c"i = / r'^drdn tr 



I A_ (H, - 4) + A_//,$, 
i=o 



(47) 



where all the are the values obtained by diagonalizing the Ihs of the j-th order energy 
correction equation. 



IV. CONCLUSIONS 



In this letter we have developed a systematic perturbation of the Shcrodinger equation 
for a bound state of Qq at rest in l/mq by using just once the naive Foldy-Wouthuysen-Tani 
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transformation on the heavy quark sector. The effects of the FWT transformation is, as 
can be seen from Eq. (^), to treat the heavy quark as a non-relativistic color source. It is 
shown that all the positive components of wave functions at each order of perturbation are 
proportional to the lowest order solution, $^ except that cf^ = 0. There also appear negative 
components of wave functions for heavy quark sector, which can also be perturbatively 
solved. As one can easily notices, the lowest order equation, Eq. (^), can be derived by 
any method. That is, we do not need to introduce the FWT or any other transformation 
to obtain this equation. The higher order terms in I/itlq, however, reflect what kind of 
transformation is introduced and we believe that any transformation which utilizes the 
properties of a heavy quark will give a better perturbation as has been shown by HQET 
than that without a transformation. 

To make formulae simple in this letter, all the gamma matrices of light and heavy quarks 
are multiplied from left with the wave function, To make the wave function a true 
bispinor, the Schrodinger equation should be transformed by a charge conjugation operator 
Uc = ilqlq for a heavy quark. That is, Ucip^ becomes a bispinor, UcHp^tU^ is a hamil- 
tonian of interest, and gamma matrices of a light anti-quark is multiplied from left while 
those of a heavy quark from right. 

One example for model dependent interaction terms is, e.g., given byllB 

PgpQ S{r) + |l - i [a, ■ + (a, ■ n){aQ ■ n)]j V{r), (48) 

where n = f/r. This model is being studied by using our formulation, in which masses and 
eigenfunctions of D, D*, B and B* are calculated including all those corrections mentioned 



in Sec. Ill and spin-flavor symmetry at the lowest order is realized by a special quantum 



number related to the heavy quark spin.i Our formulation gives the same results as those 
ini if negative components are neglected and if only the terms up to 1/ rnq in the hamiltonian 
are taken into account. 

In order to calculate, e.g., the semileptonic weak decay B ^ D + iu, one needs to boost 
the D meson even if the B meson is at rest. The formulation obtained in this letter does 
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not give any prescription how to construct a boosted wave function. We will propose some 
method in a future paper to compute the form factors by using wave functions calculated 
by the method proposed in this letter. 



ACKNOWLEDGMENTS 

The author would like to thank K. Akama and T. Morii for helpful discussions and also 
the theory group at Institute for Nuclear Study for a warm hospitality where a part of this 
work was done. 



13 



REFERENCES 



^M. Voloshin and M. Shifman, Sov. J. Nucl. Phys. 45, 292 (1987) and 47, 511 (1988); 
N. Isgur and M. Wise, Phys. Lett. B232, 113 (1989) and B237, 527 (1990); 
E. Eichten and B. Hill, Phys. Lett. B234, 511 (1990); 
B. Grinstein, Nucl. Phys. B339, 253 (1990); 
H. Georgi, Phys. Lett. B240, 447 (1990); 

A. Falk, H. Georgi, B. Grinstein and M. Wise, Nucl. Phys. B343, 1 (1990). 

2M. E. Luke, Phys. Lett. B252, 447 (1990); 
J. G. Korner and G. Thompson, Phys. Lett. B264, 185 (1991); 
T. Mannel, W. Roberts and Z. Ryzak, Nucl. Phys. 368, 204 (1992); 
S. Balk, J. G. Korner and D. Pirjol, Nucl. Phys. B428, 499 (1994); 

^ See for review, 

N. Isgur and M. B. Wise, B Decays, (ed. S. Stone, World Scientific, Singapore, 1992); 
H. Georgi, Perspective in the Standard Model, (ed. R. K. Ellis, C. T. Hill and J. D. Lykken, 
World Scientific, Singapore, 1992). 

^L. L. Foldy and S. A. Wouthuysen, Phys. Rev. 78, 29 (1950); 
S. Tani, Prof. Theor. Phys. 6, 267 (1951). 

^ J. Morishita, M. Kawaguchi and T. Morn, Phys. Rev. 37, 159 (1988). 

^ A talk given by T. Matsuki at ConGnement'95 and also at Wein'95 held at RCNP (Osaka 
Univ.), TKU-95-1 preprint, Sept. 1995, to be pubhshed in the Proceedings of Wein'95; 
T. Matsuki and T. Morii, TKU-95-3 preprint in preparation. 

■^S. N. Mukherjee, R. Nag, S. Sanyal, T. Morn, J. Morishita and M. Tsuge, Phys. Rep. 231, 
201 (1993). 



14 



